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Abstract 

We study the relationship between entanglement and spectral gap for local Hamiltonians in one di- 
mension. The area law for a one-dimensional system states that for the ground state, the entanglement 
of any interval is upper-bounded by a constant independent of the size of the interval. However, the 
possible dependence of the upper bound on the spectral gap A is not known, as the best known gen- 
eral upper bound is asymptotically much larger than the largest possible entropy of any model system 
previously constructed for small A. To help resolve this asymptotic behavior, we construct a family of 
one-dimensional local systems for which some intervals have entanglement entropy which is polynomial in 
1/A, whereas previously studied systems, such as free fermion systems or systems described by conformal 
field theory, had the entropy of all intervals bounded by a constant times log(l/A). 

1 Introduction 

In many local quantum systems, the long-distance behavior of correlation functions is a good diagnostic 
of the critical properties of a system near a phase transition as the excitation gap A tends to zero. 
Entanglement provides a different diagnostic, which gives complementary information about critical 
behavior, and is important for our ability to simulate these systems on classical computers, as described 
later. 

To quantify entanglement, we define the entanglement entropy of a region A with the complementary 
region B (assuming the global state on A U B is pure) as the von Neumann entropy S(pa), with pa 
the density matrix of region A (i.e., the global state traced over region B). A state of a spin system 
that obeys an "area law" has the property that the entropy of a region (in one or more dimensions) is 
bounded by the surface area of the region instead of the usual maximum of the volume of the region. 
(See PP for a review of area laws and references to a full cross-section of work in that area.) Area laws 
are conjectured to hold for the ground state of any gapped Hamiltonian spin system with short-range 
interactions in finite dimensions, but have not been proved in general. For one-dimensional systems, the 
area law simply states that the entropy of an interval along the line is bounded by a constant, no matter 
how large the interval. 

At a critical point, it is possible that the entanglement, like the correlation length, can diverge. 
However, if we are merely near a critical point (either because of additional Hamiltonian terms or because 
we are considering a finite-size system), the gap A becomes finite, and in one dimension, an area law 
must hold [2]. The entropy can thus provide another way to identify a critical point. Note that diverging 
entanglement and diverging correlation length do not necessarily occur together; on the one hand, there 
exist states with all correlation functions short-range but with arbitrarily large entanglement 3,4, while 
on the other hand there are critical points in two dimensions where the correlation functions become 
long-range but the area law is still obeyed [5]. 
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In one dimension, the entanglement entropy has important implications for our ability to simulate 
quantum systems on a classical computer. Intuitively, so long as the entanglement entropy remains 
small, the system is "classical" in some sense. More precisely, we know that whenever the system obeys 
an area law for the Renyi entropy, S a = (1 — m (tr(p5i)), f° r some a < 1 (which implies an area 
law for the von Neumann entropy), there exists a matrix product state approximation to the ground 
state with an error inverse polynomial in the bond dimension [6]. This approximation lies behind the 
powerful density-matrix renormalization group algorithm [7] , as well as behind Vidal's algorithm for time 
evolution [B]. 

Hastings [2] proved that a general one-dimensional gapped Hamiltonian spin system satisfies an area 
law. However, the bound is probably not tight: 

^(region) < S max = c £' In f ' In D 2 e ' ln D , (1) 

where Co is a constant of order 1, D is the Hilbert space dimension of a single spin in the system, and 

£' = 6max(2«/A,£ c ), (2) 

where v and £ c are parameters from the Lieb- Robinson bound [91 1101 ITT1 112] . and A is the spectral gap 
of the Hamiltonian. The quantity v is the Lieb-Robinson velocity, and is roughly equal to the size of 
the energy scale of the largest individual terms in the Hamiltonian times the interaction range. The 
Lieb-Robinson velocity sets an upper bound for the group velocity of excitations, so that for a local 
operator O, the operator 0(t) — exp(iHt)0 exp(— iHt) can be approximately described by an operator 
supported on the set of sites within distance vt of the support of O. The quantity £ c is a length of order 
the interaction range which determines the distance which excitations can "leak" outside the light cone 
of the Lieb-Robinson bound. In particular, all of these constants can be taken to be order 1, except for 
A, which could be very small even with all other parameters constant in size; for example, by tuning a 
system through a quantum critical point the quantity A can be made small. Thus, in the bound (Q]), the 
entropy is limited by exp[0(l/A)]. 

In contrast, physical intuition suggests that entanglement should be able to reach a distance roughly 
equal to the correlation length, i.e., about 1/A. However, most one-dimensional systems considered by 
physicists in fact have an even smaller amount of entanglement. Near a conformal critical point, the bound 
in the area law diverges as only log(l/A) |13l I14j . Indeed, as far as we aware, for all one-dimensional 
systems where the area law has previously been calculated (including those near non-conformal critical 
points), the entropy of an interval is also bounded by about log(l/A), which is below the bound in eq. ((1} 
by not one but two exponentials! 

Therefore, we raise the following as an open question: given a one-dimensional system with interaction 
range and interaction strength both 0(1), how big can the entanglement be for a given spectral gap? 
Someone familiar with the previous work on the subject might well conjecture a bound of log(l/A). In 
this paper we falsify that conjecture and get closer to the naive intuitive bound of 1/A by demonstrating 
a family of systems where the entropy of some regions is poly(l/A). In particular, we construct systems 
on a pair of rings coupled along a junction, each with N 3-state spins. In the ground state, the two 
rings are highly entangled with each other, and the entropy for one ring is roughly N. The spectral gap 
A is roughly l/7V 4 log7V, so these systems have entropy for some regions equal to (— A log A) _1//4 . The 
Hamiltonians we present depend on N in two ways: First, they are associated to the ring structure, and 
the size of the rings changes with N, and second, we have some terms in the Hamiltonian which are 
polynomially small in TV. Since there are also Hamiltonian terms of size 0(1), the presence of some small 
terms in the Hamiltonian does not affect the bound (JTJ. In sec. [3] we discuss how to let the number of 
the particles in the system be much larger than N, so we can consider TV to be a separate parameter 
from the size of the system. 

We have not succeeded in obtaining the naive estimate of 1/A and we are very far from exp[0(l/A)], 
so at present neither of the authors has any firm intuition for a tight upper bound. Given progress in recent 
years in understanding the relation between correlations and spectral gap [10], hopefully this question 
will be answered too. Intuitively, an area law is fairly simple: in a gapped system, the entanglement 
spreads only over a finite length and therefore only the spins near the boundary of A are entangled with 
those outside of A. However, in general any spin within A will be entangled with other spins within A 
as well as spins in the complement of A. Thus, to turn this intuition into practice seems quite tricky as 
it leads us into the difficult question of multi-party entanglement. 
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The intuition for the construction used in this paper is somewhat similar to that used for proofs of 
QMA-completeness and universality of adiabatic quantum computation 15 , although the details are 
rather different. In particular, we imagine a computational process that generates entanglement among 
the spins and write down a Hamiltonian whose ground state is a superposition of time steps of the 
process. 

To be more precise, we imagine the spins have the three states |0), and \x). |0) and |1) are 
two states of a qubit, and most of the spins in each ring will be in this qubit subspace, but one spin 
in each ring will be in the \x) state. We can imagine the |a;} as a "hole" in a ring otherwise filled with 
spin- 1/2 particles. The rings will meet in two sites, which we number 1 and N on each ring, with sites 
2, . . . , N — 1 on each ring being connected only to points in order on the same ring. When there are 
two qubit states at the two number-1 sites, we project onto a maximally entangled state such as the 
singlet |0}l|1).r — |l)z,|0)ij. We then move each \x) state around its ring, shifting the qubit states by 
one place, allowing us to entangle another pair. The steady-state behavior of this system is thus N — 1 
maximally entangled pairs, aligned roughly between corresponding sites on the two rings, perhaps offset 
by one place, depending on the locations of the |a;} states. 

We then write down a nearest-neighbor Hamiltonian whose ground state is the superposition of such 
states over all possible locations of the two \x) states. In the ground state, the entropy of one ring 
taken alone will be slightly over N — 1, since we broken all N — 1 EPR pairs, but we will show that the 
spectral gap of the Hamiltonian is A = Q(l/N 4 log iV). Thus, the entropy of just the left ring, say, is 
fi((-Alog A)~ 1/4 ), a polynomial in 1/A. 

The ground state differs from those used in the previous QMA-completeness constructions in that 
there is an \x) state in each ring, and we take the superposition over possible positions of the pair, 
whereas in the previous constructions, the ground state was a superposition over time slices. Essentially, 
the previous Hamiltonians were quantum walks on a line, whereas ours is a walk on a two-dimensional 
graph. 

While the Hamiltonian we write requires fine-tuning, it is not too unnatural from the physics point of 
view. Regarding the state |0) as representing an electron with spin up, the state |1) as an electron with 
spin down, and \x) as a hole, within the rings our Hamiltonian is simply the t — J model Hamiltonian 
of condensed matter physics with t / and J = 0. Only where the rings meet do we have a more 
complicated interaction. 

Sandy Irani has shown a similar result independently [16]. The two constructions have slightly dif- 
ferent motivations and each has different advantages. We were primarily interested in the dependence of 
entanglement entropy on gap, and we have a more rapid divergence of entanglement entropy with inverse 
gap and use fewer states per site. Irani's primary motivation was to study the effect of translational 
invariance on entanglement entropy. Her construction shows that it is possible to achieve a polyno- 
mial dependence of entanglement entropy on gap with a translation-invariant Hamiltonian, whereas our 
construction explicitly breaks translation invariance. 

2 From Two Rings To An Infinite Line 

As mentioned above, our main concern in this paper is to construct a local Hamiltonian H for a system 
of 2N 3-state spins. It will consist of two rings containing N sites each which are connected at two 
adjacent points, and we will show that for the ground state of H, the entanglement of the left ring with 
the right ring is 0((— A log A) _1//4 ) (where A is the spectral gap of H). It may not be immediately 
obvious how this relates to the area law of an infinite one-dimensional system, so before we discuss the 
technical details of the construction, let us first explain how the result immediately carries over to a 
system of spins on a line with the property that the entanglement of some regions is polynomial in the 
inverse gap. 

First, note that we can squash the two rings down into a one-dimensional "brick" composed of a total 
of N 9-state spins, as illustrated in figure[T] Note that 3-state spins which are neighboring in the rings get 
mapped to neighboring 9-state spins in the brick, so the resulting Hamiltonian is still a nearest-neighbor 
Hamiltonian. While at the place where the two rings meet, the ring Hamiltonian has a four-spin term 
(see the following section), all four of those spins are contained in just two 9-state spins, so the brick 
Hamiltonian contains only two-spin terms. Note that it is important that we squash the rings this way 
and not try to connect the rings along their length, as we are interested in studying the area law based on 
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Figure 1: a) The two rings. Each black dot represents a 3-state spin, b) The squashed rings. Now, vertical 
pairs of black dots represent a single 9-state spin. 



Figure 2: Part of an infinite system of bricks. 

the entropy of a single ring by itself, and we therefore want the surface area of a single ring to be small. 
In the brick, the surface area of what used to be a ring is simply 2 spins: the 9-state spin consisting of 
the original sites 1 and N, and the 9-state spin consisting of sites N/2 and N/2 + 1. If we had joined the 
rings length-wise, the surface area would be N, the size of a ring, which would be no good. 

Now we have a finite truly one-dimensional system whose ground state has the property that the 
left half of the brick has high entanglement with the right half. We can take an infinite line of 9-state 
particles and chop it up into identical bricks of this sort (figure [5]). If we take a Hamiltonian which 
is simply the sum of the brick Hamiltonians on each piece, with no terms interacting distinct pieces, 
the ground state of the resulting Hamiltonian is simply the tensor product of the ground states of the 
individual bricks. There are many excited states, but the lowest-lying excited states are superpositions 
of states for which all but one brick is in the ground state, and the remaining brick is in its first excited 
state. The gap of the line Hamiltonian is thus just equal to A, the gap of the original ring Hamiltonian 
(i.e., fi(l/JV 4 logiV)). Let us examine an interval A which begins at the center of one brick and ends at 
the center of another brick. Bricks which are completely contained in A do not contribute to the entropy, 
but the two bricks which are cut in half each contribute about N units of entropy, so the entropy of A 
in the ground state is roughly 2N. In particular, the entropy of A is fl((A log A) -1 / 4 ). 

Of course, some regions have very low entropy. Indeed, an interval containing only complete bricks will 
have entropy 0. This does not affect the area law for the system, which simply provides an upper bound 
on the entanglement of any region. Some regions have high entropy whereas others do not. Nevertheless, 
if one finds this result unsatisfying, one may always consider a larger system constructed by two adjacent 
lines of 9-state particles, each constructed as above, but offset by half a brick (as in figure [3]). Now any 
interval will cut off on each end at least one-quarter of a brick on one of the two lines. We can combine 
the two lines into a single line of 81-state particles, and the resulting Hamiltonian has the property that 
for its ground state, any interval of at least N spins will have entropy fi((— A log A) _1//4 ). 



Figure 3: Two layers of offset bricks. 
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3 The Two-Ring Hamiltonian 

We will write our Hamiltonian H as a sum of various terms: 

H = H l +H r + H b + Hp + H v . (3) 

The terms Hl and Hr cause the spins (containing a qubit state) to hop around the left and right rings, 
respectively. Hb handles the hopping of the \x) "holes" between the sites numbered 1 and ./V on each 
ring. In order to keep the two rings in sync, we require that the holes hop in unison when they hop 
between sites 1 and JV, and Hb ensures that this happens. Hp is the projection Hamiltonian, which is 
responsible for ensuring that any pair of qubits at site 1 on the left and right rings are in fact in the 
singlet state. H x is a self-energy term for the \x) hole states, ensuring that the ground state has exactly 
one | a;) in each ring. 
More specifically, let 

sf /fl = El afe >M+i<H- (4) 

a,b 

L I R 

We work modulo TV, so site N + 1 should be considered the same as site 1. S i swaps the particles in 
sites i and i + 1 for the left (L) or right (R) ring, regardless of the values of the spins at those sites. Thus 
we can define 

F L /R = S L /R]X) L /R{xl (5) 

which swaps an \x) at site i with a state at site i + 1, and annihilates any state which does not have an 
\x) at site i. 
We then let 

N-l 

H L/ a = E i\*)i /R W + K'fW K L/R {Ft'*)') • (6) 

i=l 

Hl thus produces an energy penalty for having the \x) in the iih site in the left ring, but cancels that 
energy penalty if the \x) hops to the (i + l)th site, swapping with the state there; similarly, there is a 
penalty for having the |a;) in the (i + l)th site of the left ring which is cancelled by letting it hop to the 
ith site. Hl thus tends to produce a superposition of \x) over sites 1 through TV in the left ring, with 
the qubit states shifting appropriately to make room, and Hr favors a similar superposition over sites 
in the right ring. Note, however, that Hl and Hr do not have terms allowing |a;) to hop between sites 
1 and N. 

Instead, we let Hb provide those terms: 

H B = \xx) L N R {xx\ + \xx)i R (xx\ - F^F^ - (F^ (F^ . (7) 

Hb has a structure similar to Hl and Hr, but it instead provides a penalty only if the \x) states for 
the two rings are both at site 1 or both at site TV, and then gives a cancelling energy bonus if they hop 
together to the other pair (i.e., both hop from site JV to site 1 or from site 1 to site N). In conjunction 
with Hl and Hr, this means that the \x) states in the left and right rings can wander separately in the 
two rings, but cannot complete a full circle around a ring except by hopping together. 
Hp is defined as a projector: 

Hp = (I-\x)i{x\)(I - \x) R (x\) - (8) 

which produces an energy penalty for any state which has a qubit state in site number 1 for both rings, 
unless the two qubit states form a singlet state I* - } = (|01) — |10})/%/2. 
Finally, let 

N N 

H^ R = Vi - Vi E /R (*\ + ^E \ xx )i!iU xx l ( 9 ) 

i=l i=l 

and 

H V = H^ + H R . (10) 

We choose Vi and V2 to be positive constants, so Hv will produce an energy bonus proportional to the 
total number of \x) states in the two rings, but an energy penalty if any two \x) states are adjacent in a 
single ring. Since Hl, Hr, and Hb will tend to produce superpositions of different locations for the \x) 
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states in each ring, if there is more than one jx) within a given ring, Hl, Hr, and Hp will favor states 
where some \x) states are adjacent, which will then get an energy penalty due to the V2 term. By tuning 
V\ and V2 appropriately, we will ensure that the ground state has exactly one \x) state in each ring, as 
desired. The constant term 2Vi then gives that ground state energy. 

Note that the terms Hl, Hr, Hb, and Hp are all non-negative operators. However, H v and Hy are 
not, which means we will need to exercise some caution in bounding the energy gap. 

4 The Ground State and the Gap 

We now turn to determining the structure of the ground state of H and its spectral gap. We shall prove 
the following: 

Theorem 1. Suppose V\ = 1/N 4 and V2 = 1. Then 

a) H has a nondegenerate ground state \g), 

b) The spectral gap A of H is tt(l/N 4 log N), and 

c) S(p L ) = S(p R ) >N-1, where p L/R = Tr R/L \g)(g\. 

The first thing to notice is that no term of H changes the number of \x) states in a ring. Therefore, 
we can treat the Hilbert space as a direct sum of subspaces with different numbers of \x) states in the 
two rings and analyze them separately. Let the subspace with a \x) states in the left ring and b \x) states 
in the right ring be H a ,b- We start with rli,i, which we claim contains the global ground state. We 
therefore analyze the ground state (which will have energy 0) and gap within this subspace. Then we 
show that all states in the other subspaces have energy at least 0(1/N 4 ). 

4.1 One I x) in Each Ring 

Within the subspace Hi,i, the Hamiltonian term Hy becomes uniformly 0, so we can ignore it. We 
are thus left with H = Hl + Hr + Hb + Hp, and we wish to find the ground state and gap of this 
Hamiltonian on 

First, let us define some convenient basis states. We start with states which have both \x) states in 
site N in their respective rings, and define 

|{Qi}} = \ai)i R \a.2)2 R ■ ■ ■ \a N -i)N-i\xx)N R , (11) 

where |e*i) is a Bell state, i = 1, . . . , TV — 1. The state thus has the two qubit states at the jth 

position in each ring sharing the particular entangled state \otj). 

Now define M a ,b to be an operator that shifts \x) to site a in the left ring and site 6 in the right ring: 

M a ,b = (SaSa+i ■ ■ ■ S^-l^j (^S r S r +1 ■ ■ ■ S R -ij ■ (12) 

(Make the convention that if a or b is equal to N, do not include the L or R operators, respectively. Thus 
Mn,n = I-) Then M ay t\{cti}) has the \x) states at sites a and b in the left and right rings, respectively, 
but some of the entangled states |a») are now misaligned: If j < min{a,6} then the two sites j share 
\ctj), and if j > max{a,6}, then the two sites j + 1 share \ctj). However, if a < j < b, then the state \cej) 
is shared between site j + 1 in the left ring and site j in the right ring, whereas if 6 < j < a, then \ctj) is 
shared between site j in the left ring and j + 1 in the right ring. Note that for Mn,i\{oh}), all pairs are 
misaligned by one. Thus, even though S r Mn,i\{ch}) has both \x) sites back at site N, the entangled 
pairs are misaligned, so it is not a state of the form | {/?;}} for any {Pi}. 

Let R = M\ 1 SnS§. Then R\{aa}} also has both \x) states back at site N, but since the two \x) 
states have both circled the ring once, the pairs remain in sync. They have all been shifted by one place: 

R\{ ai }) = \{fr}), (13) 

with \/3i) = |a i+ i> (i = l,...,N-2) and \/3 N -i) = \ai). 

We now look at the action of the Hamiltonian on these states. The states M ay t\{cti}) form a basis for 
Hi,i, and Hp is diagonal in this basis, with eigenvalue for any state of the form Mi^\{ai}), M a ,i\{cti}) 
(for any a, b), or for M a ,b\{cti}), with a, b > 1 and |ai) = j* - ). Hl and Hr will interact different 
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states Ma,b\{ai}) with the same list {a;} of Bell states. Hb will couple Mi,i|{ai}) with R |{a;}}, so the 
combination Hl + Hr + Hb couples all states of the form M a ^\{cti}) with the list of Bell states {oti} 
fixed up to cyclic shifts. Suppose the list {ai} has period p under cyclic shifts (i.e., R p \{ai}} = |{o!i}))- 
Then it follows that 

\S({ ai })) = -^^^{0,}) (14) 

VP N r = a.b 

is an eigenvector of Hl + Hb + Hb with eigenvalue 0, and all states annihilated by Hl + Hb. + Hb are 
of this form. We note that the Hilbert subspace Hi.i further decomposes into invariant subspaces of H , 
each of which can be labeled by a list of Bell states {on}, with two lists equivalent if they differ only by 
a cyclic shift on i. We will refer to the subspace associated to the list {on} as rl({ai}); we choose some 
convenient representative list {ai} from the set of equivalent lists to label each H({ai}). 

The state \g) — \S({ai})) with \at) = |^ - } for all i also has eigenvalue for Hb, as every term 
in it has either \x) in site 1 for at least one ring or has a singlet shared between the two 1 sites. All 
other states ^({a^})} do not have this property: Even if |ai) = l^ - }, some \ai) 7^ |^~), and therefore 
HpR 1 ~ 1 \{ai}} = 1. \g) is thus the ground state within (and, we will show, within the whole Hilbert 
space), whereas we claim any other state has substantially higher energy. 

Theorem 2. If \ip) 6 H({ai}) with \on) 7^ I 1 ! 7 } for at least one value of i, then 

(4>\Hl + Hr + Hb + H p \^) = fi(l/p 2 7V 2 log/V), (15) 

where p is the periodicity of {at} under cyclic permutations. 

Since p < N, the theorem implies that the gap of the full Hamiltonian within Hi.i is Q(l/N 4 log TV). 

As a warmup to proving this, we first consider Hl + Hr + Hb and prove that it has a polynomial 
gap. While we will not use precisely this proposition in the proof of Theorem the proof uses a very 
similar technique, and the main lemma for Prop. Q]will also be central to proving Theorem [21 

The Hamiltonian Hl + Hr + Hb can be thought of as a quantum walk on a graph. In particular, let 
us restrict to the invariant subspace H({cti}), where {ai} has period p under cyclic shifts. We then get a 
graph G p of pN 2 nodes. The states |o, 6, r) = M a ^ b R r \{ai}} can be considered as the nodes of the graph, 
which consists of p N x N grids, with the (N, N) site of one grid connected to the (1, 1) site of the next 
one, and the (TV, TV) site of the pth grid connected to the (1, 1) site of the 1st grid. (See figure[3]) Edges 
of the graph correspond to states which are coupled by a single off-diagonal term of the Hamiltonian. 

Proposition 1. The gap of H L + H R + H B within %({«»}) is Q.(l/p 2 N 2 log TV). 

Proof of proposition. Hl is a simple quantum walk on a line of length N. It can be explicitly diagonalized, 
with eigenstates 

Wk,b,r)=J]- J £cos(^(a-l/2)) \a,b,r) (16) 
0=1 ^ ' 

A, =2(1- cos ^) (17) 

for k = 1, . . . , N — 1. fc = Ois also possible (although with a different normalization constant), and gives 
the ground state \tp , b, r) — (\/\J~N) "}2 a \a, b, r) , with energy 0. Similarly for Hr. Since a single grid has 
the two coordinates completely uncoupled, the eigenstates \ipu,r) of Hl + Hr are similarly simple, with 
eigenvalues Afc; = Afc + A;. The ground state of Hl + Hr (for fixed r) is \ipoo, r) = (1/N) ~}2 a b \a, b, r). 
For general k,l, let us write \ipkl,r) = J2 a ,b fki( a , b)\a, b, r). Then 

fki{a, b) = C kl cos fea - 1/2)) cos ^(6 - 1/2)) , (18) 

where Cm is a normalization constant: Cm = 1/N for (k,l) = (0,0), Cm = V2/N if exactly one of k, I 
is 0, and Cm = 2/7V otherwise. 

The ground state of Hl + Hr + Hb is again simple: 

j p-l j p-l N 

Iground) = -=^2\1>ao,r) = — — ^ ^ |a,6,r>. (19) 

V r r = V r r = o a ,6=l 



and eigenvalues 
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\N,N,2) 



11,1,2), 



\N,N,1) 



\1,N,0) 



11,1,1), 



\N,N,0) 



/ 1,1,0) JV.1,0) 



Figure 4: The graph of basis states within H({ai}). 



However, the higher eigenstates are more complicated. We can simplify slightly by noting that the 
graph has translation invariance under the action of R, so the eigenstates of Hl + Hr + Hb will also 
be eigenstates of R. R has eigenvalues exp(2-7rim/p) (since R p = I on the subspace H({ca})). Any 
eigenstate of J? can be written as 



p-l N-l 
V P r=0 fe,I=0 



(20) 



with m = 0, . . . ,p — 1. Let us consider the expected energy (V>|-ffL + JTr + Hb\iP) of any state \ip) which 
is orthogonal to the ground state | ground) and is an eigenstate of R. If m — 0, coo = (because the 
state is orthogonal to the ground state), so 



Oj\H l + H r \^j) > Ai = 6(l/iV 2 ). 



(21) 



Therefore, we should search for the first excited state among states with nonzero m. In particular, we 
will show that for any state with m > 0, either the state is dominated by the ground state of Hl + Hr 
(i.e., coo is very near 1), in which case it does not satisfy Hb very well, or there is a significant mixture 
of higher eigenstates of Hl + Hr. 

For the proof of Theorem we will need a stronger version of the result using unnormalized states. 
The necessary result is contained in the following lemma: 



Lemma 1. Let 



Then 



where m! = min(m,p — m) 



p-l AT-l 



" mr/p \i>ki,r). 
\H L +H R + H B \^) = Q(m' 2 \c ao \ 2 /p 2 N 2 log AT), 



(22) 
(23) 



When \ip) is normalized, it is easy to see that {iP\Hl + Hr\iP) is large (Q(l/N 2 )) unless |coo| is order 
1, so the proposition follows immediately from the lemma. 

□ 
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Proof of lemma\J] Let us consider the expectation value of Hb for m > 0: 

p-l N-l 



(#ff S |V>> = -£ E Cwc^e^-^^^.r'IffB^w.r). (24) 



P 

r,r'=0 k,l,k',l'=0 



Since 



(<l> kn ,,r'\H B \<p kl ,r) = 6 ry [/w(l, l)/^;/(l, 1) + fki{N,N)fk'i'(N, N)] 

- 5 r y^f k i(N,N)f^,(l, 1) - <V +1 / fei (l, l)fr v (N,N), (25) 

we have, focusing on the (k, I) = (0, 0) terms, 

(Vl^slV) = I coo 1 2 [2 - 2cos(27rim/p)] + (tp'\H B \^'} 



N 2 
+ 



^ E c^cSo [/ fei (l,l)(l-e 2 ™ /p ) + / fei (iV,iV)(l- e - 2 ™/ p )] 
(fe,0#(o,o) 

+ i E coo^f [AV(l,l)(l-e- 2,rim/j ') + /^ i 'W^)(l-e 2,rim/p )], (26) 



(fc', !')/(<>,()) 



where 



l^> = -^E E c fe! e 2 " m ^|V> W! r). (27) 

V P r = (fc,i)^(0,0) 

Since i/s is positive semidefinite, (iP'\Hb\iP'} > 0, and 

(i>\H B \ij}> ^jc t ,oi 2 [2-2cos(27rm/p)] (28) 



+ 2Re |^f E c fei cS [/ fc; (l,l)(l-e 2 ^) + / fc; (iV,iV)(l- e - 



1 



Inimj'p \ 



= ^|coo| 2 [l-cos(2^m/p)] (29) 
+ Jf E Re(c fci c 00 )C fci cos(7rfc/27V) cos(irf/2JV)[l - cos(27rm/p)][l + (-l) fc+i ] 

fc,i^(0,0) 

+ E Im(c W cSo)CAi cos(7rfc/2iV) cos(7rZ/2iV) sin(27rm/p)[l - (-l) fc+ '] 

fc,i^(0,0) 

> I 1 - cos(27rm/p)] (30) 

g ^ g ^ 

_ Jp E \ c ki\\coo\[l - cos(2-Kin/p)] - E |cfci||coo||sin(27rm/p)| 

fc-f-Z>0 even k-\-l odd 

= ^i[l-<-..M-.'^ /-) ( |,„„| - I V |o.,|- l|<-<.K-m- / ,i| ( ; > I ) 



[1 - cos(27rm/p)] I jc 00 | - 4 E l c «l _ 4 I cot(7rm/p)| E 

\ fc + Z>0 even k + l od 



fc + Z>0 even fc + Z odd 

For < m/p < 1/2, | cot(-7rm/p)| < p/(mn), so 



(^IV*) > ^[l-C0B(2irm/p)] 



|coo|-4 E E M 

fc + Z>0 even fe + Z odd 



(32) 



For 1/2 < m/p < 1, we can just consider m! = p — m, in which case | cot(-7rm/p)| = | cot(-7rm'/p)| an d 
cos(27rm/p) = cos(27rm'/p), and we can just substitute m for m in equation ((32}. Otherwise, we let 
m! = m. 

Therefore, one of two things must be true. The first possibility is that the sum is dominated by coo: 
cool - 4 E M-4-7- E M>|coo|/2, (33) 

* — ' 771 TV A — J 

fc + Z>0 even fc + Z odd 
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m which case 

m |c o| 



l^slV) > [i - cos(W/p)] = n (- 



The second possibility is that there is a non-negligible contribution from the higher k, I eigenvalues: 

V 



(34) 



cool -4 V |c fci |-4-^- V |cfci| < |c o|/2, (35) 



m'7r 

k + l>0 even fc + Z odd 



which implies that 



However, 



Note that 



for < k < N, so 



\ckl\ =n(m'|coo|/p). (36) 

M#(0,0) 

JV-1 

(V»|fT £ + H R \1>) = J] ( Afc + X ')\ c "i\ 2 - ( 37 ) 

fc,Z=0 

A fc =2(l-cos(^))>^ (38) 
MHl + HrMZ-^ ik 2 +f)\c kl \ 2 . (39) 

fc.i^(O.O) 

Under the constraint that J^ fc i7 ^ j |cfci| is fixed (as suggested by (|36|> ). the RHS of ((39j) is minimized 
when |cfc;| = A/(k 2 + I 2 ) (this minimization can be done using Lagrange multipliers, for instance), with 

E fea4]2 =«(m'|coo|/p). (40) 

fc,Z#(0,0) 

For large iV, £ l/(fc 2 + Z 2 ) = 6 (log AT), so A = Q(m'|coo|/plogiV). Thus, 

(V|fTz, + ffalV) > «(m' 2 |c o| 2 /p 2 iV 2 log 2 iV) ^ _i_ = fi( m ' 2 | Coo | 2 / p 2 iV 2 log7V). (41) 

That is, regardless of the size of |coo|, we have that {%I>\Hl + Hr + Hb\iP) = O.(m' 2 \coo\ 2 /p 2 N 2 log AT), 
proving the lemma. 

□ 

We next turn to proving Theorem [2] The proof is somewhat similar: The ground state | ground) of 
Hl + Hr + Hb is completely delocalized in the \a, b, r) basis, so in particular has substantial energy for 
Hp. In order to cancel the bad \a,b,r) sites, we will need in to add a nontrivial contribution of higher 
k, I, and m values, which will increase the energy of Hl + Hr + Hb- 

Proof of Theorem^ Let 

1 p_1 

IfeJ^E^^'^ ( 42 ) 
VP r= o 

and 

JV-l 

\?P) = Y, a klm\lpkl m ) ■ (43) 
k,l,m=0 

We start by considering |(V>fcZm|-ffp|V'ooo}|- To calculate this, we expand |V>fcZm) and |^>ooo) m the 
\a, b, r) basis, since Hp is diagonal in that basis. Then {a' , b' , r'\Hp\a, b, r) is unless a — a' , b — b' , and 
r = r' . It is also if a = 1 or 6 = 1. In all other cases, (a, b, r\Hp\a, b, r) = if \a r ) = |*& _ } and 1 
otherwise. Let 

R={r s.t. \a r ) + I* - )}. (44) 
That is, R is the set of "bad" r, that contribute to a nonzero energy for Hp. 
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Then 



\(-ipki m \Hp\ip 00 o)\ 



e 27 " mr /p {ip k i,r'\H P \ tpoo , r) 
^ J e ~ 2lTlrm/r ' '(ipki,r\Hp\ipoo,r) 



1 

<^|(Vfci,0|i?p|Voo,0>|. 



r6fi 



(45) 



(46) 



(47) 



In the last line we assume, without loss of generality, that G R. It is also worth noting that for m = 0, 
we have equality in the last line. Furthermore, 



1 

|<Voo,0|ffph/>oo,0>| = -^ (a',b',0\H P \a,b,0) 



:,b',b=l 



(N — l) 2 
N 2 



When i / 0, 



|(Voi,0|Hp|Voo,0)| = 



^2 

AT2 



f; cos (£(&' - 1/2)") (a',6',0|iJp|a,6,0) 



V2(AT - 1) 
N 2 



£ cos (£(6-1/2) 

6=9 V 



But note that we have an orthogonality relationship for cosines, so that 



£>(£(6-l/2))=0, 

6=1 V y 



KV<H,0|ffp |^oo,0)| = V ^% 1) |cos( 7 r;/2iV)| 



when I / 0. Similarly, when k ^= 0, 



|(Vfco,0|ffp |Voo,0)| = ^ ( ^ 2 ^ |cos(7rfc/2JV)|. 



Finally, when k,l =fc 0, 



|(^i,0|Hp|Voo,0)| = 



AT2 



We use another orthogonality relationship: 



°=E 



a, 6=1 
N 



Trk 

at 1 



E 



7l7 



^(a-1/2)) cos (£(6- 1/2) 



(a- 1/2)) cos (-(6- 1/2) 



- E cos (t£(— 1/2)) - (£(6- 1/2)) cos (^(a - 1/2)) cos 

a, 6=2 ^ ' ^ ' a = 2 ^ ' ^ 

+ E «» - 1/2)) cos (^) + cos (|* ) cos ( J* 

= E cos (^(a- 1/2)) cos (£(6- 1/2) 



7rfc 
2iV 



2N 



(48) 

(49) 
(50) 

(51) 
(52) 

(53) 
(54) 

(55) 



(56) 
(57) 
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where in the last line, we have used the orthogonality relationship (|51[) for each of the two single sums. 
We then get 



| (ip k i,0\Hp\ip oo , 0) = — cos ( — ) cos ( — 



2N 



2N 



(58) 



We now turn to considering (ip\H P \ip). Let = E(k,i,m)^(o,o,o) a kim\^ki m ), so \ip) = aoool^ooo) + 
ji//). Since Hp is positive definite, 



\(ip\Hp\i/)}\ >|a 00 o| 2 W'ooo|#p|V'ooo} + 2 Re{a 000 (ip'\H P \ip 000 )) 



1000 



pN 2 



+ 2Re a 000 Ofcim(' i /'feim|-ffp|V'000) 

(fc,i,m)#(0,0,0) 



>|aooo| 2 ^i)!_ 2 E 



>|aooo| 



pN 2 



\R\(N-lf 
P N 2 



| aooo 1 1 a k im 1 1 {ipkim | Hp | i/'ooo) | 



(59) 
(60) 

(61) 



(fc,!,m)^(0,0,0) 



2E 



2^2 



N 



~i E 



ClfcOn 



fc>0,m>0 



(iV-1) 2 

!>0,m>0 V y 

I/? (N — l) 2 

|aooo| rrr, [loaool — X — Y — Z], 

pl\ * 



E 



\0>klr, 



k,l>0,m>0 



where 



X 



2V2 s - 



2^2 v^i 4 
afeool + jj—j 2^ l a ° i() l + (N-l) 2 ^ |CtMo1 
fc>o ;>o ^ ' k,i 



2^2 sr^ 

N-l 2-^ 



2^2 v - 
N-l^ 



,l>0 

4 



O-Olr 



(N-l) 



1)2 E 



Oklr, 



Z =2 |aoon 

m>0 



(62) 
(63) 

(64) 

(65) 
(66) 



There are four cases: 

• If |a oo| X Y Z > 1/2: Then |a oo| > 1/2 as well, so \{^\H P \^)\ > lR} ^} f = 0(l/p). 

• If X = 0(1): At least one of the sums E fc >o l°fcoo|, -^h" E;> l«oio|, and (N ^y2 J2 k ,i>o K« 
must also be 0(1); the others may be smaller. These three sums are all averages, and at least one 
of the terms in an average must be greater than or equal to the mean value. Thus, there exists 
(Aid, Jo) / (0,0) with |a fco , o0 | = 6(1). But 



{i>\H L + H R \i>) = ]T \a klm \ 2 \ kl > \a kolo0 \ 2 n(l/N 2 ) = fi(l/A 2 ). 



(67) 



k,l,m 

• If Y = 0(1): At least one of the sums ^ J2 m>0 E fe > \ a kOm\, j^j E m>0 E i>0 \aot m \, and 
(jv^rp- E m>0 Efe,!>o \ a ki™\ must be We h ave 

t2\ 



\H L + H R \ip}> ^ ]T|a fc!m | 2 Q(l/A 2 

(k, 0^(0,0) m 



(68) 



When E l°fcim| is fixed (summed over some subset of (k,l,m)), we minimize the sum E l a w™| 2 
(over the same subset of (k,l,m)) by taking all |afc; m | to be equal (for (k,l,m) in the subset). In 
this particular case, we should choose |afe; m j = 0(l/(p — 1)) for one of the sets {fc > 0, / = 0, m>0}, 
{k = 0, 1 > 0, m > 0}, or {k > 0, 1 > 0, m > 0} and we find ^ fcim |a fcim | 2 = Q(N/p). (Or fi(A 2 /p) 
in the third case.) In any case, it certainly follows that {iP\Hl + Hn\ip) = Q(l/N 2 ). 
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• If Z = @(1): This is the most interesting case. We note that 

{il>\Ht + H R + H B \ip) = ^{ip m \H L + H R + H B \i) m ), (69) 

m 

with \ip m ) = J2 k . \ipklm}- By Lemma Q] 

{i> m \H L + H R + H B \i> m ) =n(m' 2 \a 00m \ 2 /p 2 N 2 \ogN), (70) 

where m! = min(m,p — m). Therefore, we wish to minimize ~}2 m i~n' 2 \aoom\ 2 under the constraint 
that E m >o l a oom| = ©(1)- The minimum is achieved (as can be shown by Lagrange multipliers, for 
instance) when |aoom| = B/m' . Since 

Em 1 / m ' 2 ^ 2 C( 2 ) = tt 2 /3, it follows that B = 6(1) and 

E m ™' 2 K)m| 2 =0(1). Thus, 

{i>\H L + H R + H B \i>) = Q(l / P 2 N 2 log N). (71) 

As it happens, case 2 cannot actually occur in the large N limit, given the normalization constraint 
that ^ \ a kim\ 2 = 1, but this is not very important since the gap is set by the fourth case. 

□ 



4.2 Zero or Multiple \x) States 

We now turn attention to the subspaces H a ,b, with either a or 6 different than 1. We can treat Hy and 
Hy separately, and show that Hl + Hy is bounded below when a ^ 1, as is Hr + Hy when 6^1. Since 
the remaining Hamiltonian terms H B and Hp are non-negative, bounding Hl + Hy and Hr + Hy is 
adequate to give a lower bound on H. Since the two cases are identical, we focus on the left ring. 

If a = 0, Hy has all eigenvalues equal to Vi and Hl is 0. If a > 1, there are some states in H a ,b 
which have a negative eigenvalue for H v . However, those states will have a positive eigenvalue for Hl, 
and we must choose Vi and V2 to give such states a positive overall energy. 

Let us restrict attention to a particular value of a; then H v + (a — l)Vi is a non-negative operator, 
which allows us to apply the following lemma (see [T7] for a proof): 

Lemma 2. Let A\ and A2 be nonnegative linear operators with null spaces L\ and Li. Suppose that 
L\ n L2 — {0} and that no nonzero eigenvalue of Ai or A2 is less than v. Then 

A! + A 2 > usin 2 (6»/2), (72) 

where 6 — 6(Li,L2) is the angle between L\ and L2. 

For us, Ai — H v + (a — l)Vi and A2 = Hl- The null space L\ of A\ consists of superpositions of 
states with a \x) states in the left ring with no two of those \x) states are adjacent. The smallest nonzero 
eigenvalue of A\ is at least Vs. (It could potentially be higher than V2, if a is large enough so that all 
configurations contain multiple adjacent \xx) pairs in the left ring.) The null space L 2 of A 2 consists of 
uniform superpositions over all possible configurations of a |a;) sites on the left ring, with some state for 
the remaining qubit sites, shuffled around appropriately for different \x) configurations. 

To find the lowest nonzero eigenvalue of A 2 , we can recognize Hl as the spin- 1/2 ferromagnetic 
Heisenberg model on a chain of length N with open boundary conditions by mapping \x) to spin up and 
a qubit state to spin down. In this case, we now consider the \x) states to represent particles hopping 
around on the chain. Unlike the usual ferromagnetic Heisenberg model, the background states with no 
particles have some structure, since the background states are actually qubits. However, the qubits are 
inert and cannot hop or otherwise interact by themselves; they can only switch places with a \x) particle 
as it hops. Therefore, the Hilbert space breaks up into invariant subspaces within which the state of the 
background sites is determined completely by the locations of the \x) particle states. Within a single 
such subspace, we have exactly the usual spin-1/2 ferromagnetic Heisenberg model, which is known to 
have gap 1 — cos(-7r/iV) = Q(l/N 2 ) [18] . (This gap considers all possible numbers of \x) states; in fact, 
the minimum gap is achieved when there is just one \x).) 

The angle 8 between L\ and L2 can be determined by simply noting that there are ( ) total possible 
|a;) configurations, but at most 

N(N-2)(N-4)(N-6)---(N-2a + 2) < fof\ ^ _ ^ 
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configurations with no adjacent \xx) pairs. This formula comes from choosing the location of a \x) sites, 
with each one chosen from possibilities excluding the sites of previously chosen \x) sites as well as the next 
higher-numbered site to any previously chosen site; then we divide by a! to eliminate the redundancy 
due to choosing the sites in order. This of course is an overestimate, as we must also exclude the next 
lower-numbered site to any previously chosen site, but then the counting becomes complicated, as the 
site below one previously chosen site can be the same as the site above a previously chosen site. Since 
an element of L2 consists of a superposition over all (^) configurations and an element of L\ can only 
contain configurations with no adjacent pairs, it follows that sin 2 8 > N 1 _ 1 . 

Putting these numbers into lemma[2] we find that on the subspace with a > 1, Hl+H v + ((i — 1)Vi has 
all eigenvalues Q.(l/N 3 ). If we choose Vi = 0(l/iV 4 ) and V2 = 1, then Hl + H v also has all eigenvalues 
0.(1/ N 3 ) on the subspace with a > 1, while on the sector with a = 0, Hl will have eigenvalue Q(l/N 4 ). 
Putting everything together, we find that the ground state is unique, and is the desired ground state \g) 
(with energy 0), while the gap in the whole Hilbert space is n(l/7V 4 logiV). 

4.3 Entanglement of the Ground State 

The ground state of H is thus 

JV JV-1 

\9) = Jf E M a , b \xx)lf <g) I*")?*. (74) 

a,b=l i—1 

as noted before. Suppose we make a local measurement to each ring that measures the location of the 
|a;) state but nothing else about the qubit states. The outcomes are uniformly distributed among all 
possible pairs (a, i>), and given a specific outcome, the residual state is 

JV-1 

A/ a , b |xx)^ B (g) \<t>-}f R . (75) 

i=l 

M at b is just a tensor product of a local unitary operator M a on the left ring with a local unitary operator 
Mj, on the right ring. By performing M\ and M^, we can therefore convert the state to IVI/ - }®^ -1 ). 
Since purely local operations convert \g) into TV — 1 EPR pairs, the original state has at least N — 1 
ebits of entanglement. That is, S(pl) = S(pn) > N — 1. Since one ring has only N 3-state particles, 
the maximum possible entropy would be iVlog 2 3, so we find that ^(one ring) = Q(N). If we rephrase 
this in terms of the energy gap, A = Q,(l/N 4 log N), we see that 5*(one ring) = fi((— A log A) -1 / 4 ), as 
claimed. 



5 Discussion 

The explicit construction in this paper shows that it is possible to have the entanglement entropy scale 
polynomially with the inverse gap. On the other hand, it was proven previously that the entanglement 
entropy cannot scale faster than exponentially in the inverse gap. A tight bound remains open. 

One route to a tight bound might be to improve the Lieb-Robinson bound. This bound is a key 
ingredient in proving upper bounds on entanglement entropy, but our construction suggests that the 
Lieb-Robinson bound does not fully capture the locality properties of Hamiltonian dynamics. Let us 
explain by an example: consider a system of free fermions on an infinite line in one dimension, with 
Hamiltonian J]V \^t +1 '4 r i + h.c. Consider the operator equations of motion for the operator ^J. The Lieb- 
Robinson bound says that exp[iff£]\l/J exp[— iHt] can be well-approximated by an operator supported 
just on the set of sites within distance VLRt of 0, where vlr is the Lieb-Robinson velocity. However, at 
least for this system, we can make a sharper statement — there is no good approximation by an operator 
supported on too few sites. The exact result for this operator is: 

exp[iHt]^l exp[~iHt] = S~] A(x, t)^t, (76) 

where 

r 

A(x,t) = / — exp[2icos(kt)]exp(ikx). (77) 
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One can check from the asymptotics of the Bessel function that for x » 2t the amplitude A(x,t) is 
exponentially small, in agreement with the Lieb- Robinson bound. However, the amplitude A(x,t) is 
symmetric under a change in sign from x — > —x. That is, the operator creates both left and right 
moving excitations. Further, while the amplitude is maximum at x close to ±2t, there is a non- vanishing 
spread in the width of |A(:r, t)| 2 : the particle is not perfectly localized at x = ±2i, but instead has a 
spread in its position, indicating an uncertainty in the initial velocity of the particle. It is possible to 
create a wavepacket that moves to the right, but in order to do this we need an operator that creates the 
particle in a superposition of different sites: Yli a (*)^l f° r some function a(i). If a(i) has finite support, 
then there will always be some spread in the velocity of the wavepacket, and in fact there will be some left- 
moving component to the wavepacket. Thus, it will not be possible to find, at least in this model, a choice 
of aii) which has finite support for which the time evolved operator exp[iHt] a{i)^\ j exp[— iHt] is 
approximately localized on a set of sites with a bounded, time-independent diameter. The Lieb-Robinson 
bound does not capture this effect, but it seems to be a real effect for every Hamiltonian that we can 
think of: operators with finite support create both left- and right-moving excitations. 

A similar phenomenon afflicts our construction. We had to go to some effort to keep the two rings 
synchronized because of the spread in the velocity of the \x), which, if left unchecked, could wash out the 
entanglement between the two rings. To do this, we only allowed the \x)s to hop from site 1 to site N in 
pairs. The fact that the term Hb only has effect when both \x)s are in the correct position is one of the 
reasons that the gap is so small. However, if we had a way of just creating an \x) excitation that moved 
to the right around the ring at a constant velocity equal to vlr, it might be easier to synchronize the 
rings, and thus to increase the gap for the same entanglement. Hopefully, then, this work will lead to a 
better Lieb-Robinson bound that captures these effects, and this in turn will lead to stronger bounds on 
entanglement in Hamiltonian systems. 
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